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Abstract. The uniform recursive tree (URT) is one of the most important models and has been success- 
fully applied to many fields. Here we study exactly the topological characteristics and spectral properties 
of the Laplacian matrix of a deterministic uniform recursive tree, which is a deterministic version of URT. 
Firstly, from the perspective of complex networks, we determine the main structural characteristics of the 
deterministic tree. The obtained vigorous results show that the network has an exponential degree distri- 
bution, small average path length, power-law distribution of node betweenness, and positive degree-degree 
correlations. Then we determine the complete Laplacian spectra (eigenvalues) and their corresponding 
eigenvectors of the considered graph. Interestingly, all the Laplacian eigenvalues are distinct. 
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1 Introduction 

The flexibility and generality in the description of natural 
and social systems have made complex networks become 
an area of tremendous recent interest [UdlGIH]. Among 
many interesting aspects, topological characterization is 
very significant for the study in network field. In the past 
ten years, there has been a considerable interest in char- 
acterizing and understanding the topological properties 
of networked systems [5] . A lot of network measurements 
have been proposed, among which degree distribution, av- 
erage path length (APL), betweenness, and degree corre- 
lations have been extensively studied, since they have pro- 
found effects on the dynamical processes taking place on 
networks, such as robustness [6j[TlISj[9] , epidemic spread- 
ing [TUIITT] , synchronization [ll^[T3"l[L4] , and games [TSlfTB] . 

The above mentioned topological characteristics focus 
on direct measurements of structural properties of net- 
works. Apart from these investigations there exists a vast 
literature related to (Laplacian) spectrum of complex net- 
works [17,18,19,20,21], which provides useful insight into 
the relevant structural properties of graphs. In fact, topo- 
logical features capture the static structural properties of 
complex networks, while spectrum provides global mea- 
sures of the network properties [17] • In the past years, 
graph spectrum has found many important applications 
in physics and other fields [3,4|. For example, the ratio of 
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the maximum eigenvalue to the smallest nonzero one of 
Laplacian matrix determines the synchronizability of the 
network [L2lll3[|T4] . On the other hand, the eigenvectors 
of Laplacian matrix have also been successfully used to 
detect community structure of networks [55]. While the 
Laplacian eigenvalues and eigenvectors have high influ- 
ence on the structural properties of networks and dynam- 
ics running on them, until now, most analysis of Laplacian 
spectra and eigenvectors has been confined to approximate 
or numerical methods, the latter of which is prohibitively 
time and memory consuming for large networks [17] . 

On the other hand, in order to describe real systems 
and study their structural properties, a wide variety of 
network models have been presented [UEIEH!] j among which 
the uniform recursive tree (URT) is perhaps one of the 
most widely studied models [23] . It is now established that 
the URT is one of the two principal models [24]l25] of a 
random graph (the second one is the famous Erdos-Renyi 
model [26]). The URT is perhaps the simplest tree and 
is built in the following way: at each time step, we at- 
tach each new node to an existing node which is chosen 
uniformly at random. It has found applications in several 
areas. For example, it has been suggested as models for 
the spread of epidemics [27] , the family trees of preserved 
copies of ancient or medieval texts [28], chain letter and 
pyramid schemes [31] > to name but a few. Recently, a de- 
terministic version [32] of the URT has been proposed to 
mimic real-life systems whose number of nodes increases 
exponentially with time. This kind of deterministic models 
have drawn much attention from the scientific communi- 
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ties and have turned out to be a useful tool [33,34,35, 
36, 37, 38, 39, 40, 41 , 42, 43, 44, 45, 46, 47J . Although uniform 
recursive tree is well understood [23lf24lf25lf27U28lJ29l [30] . 
less is known about the topologies and other nature of the 
deterministic uniform recursive tree (DURT) [32 . 

In this paper, we offer a detailed analysis of the de- 
terministic uniform recursive tree (DURT) 32J from the 
viewpoint of complex networks. We first determine accu- 
rately relevant topological characteristics of DURT, such 
as degree distribution, average path length, betweenness 
distribution, and degree correlations. We then use meth- 
ods of graph theory and algebra to calculate or estimate 
the eigenvalues and eigenvectors of the Laplacian matrix. 
We present that there is a strong relationship between the 
eigenvalues and the eigenvectors of the Laplacian matrix. 



2 The deterministic uniform recursive tree 

The deterministic uniform recursive tree is one of the sim- 
plest models. It is constructed in an iterative way [32]. We 
denote the tree (network) after t steps by XJ t (t > 0). Then 
the network at step t is built as follows. For t — 0, Uq is 
an edge connecting two nodes. For t > 1, Ut is obtained 
from Ut-i- We attach a new node to each node in Ut-\- 
This iterative process is repeated, then we obtain a de- 
terministic tree with an exponential decreasing spectrum 
of degrees as shown below. The definition of the model is 
illustrated schematically in Figure [1] 

We first compute the total number of nodes N t and 
the total number of edges E t in the tree Ut- Let n v (t) and 
n e (t) denote the numbers of nodes and edges created at 
step t, respectively. Then, n v (0) = N = 2 and n e (0) = 1. 
By construction, we have n v if) — Nt-i, thus Nt = n v (t) + 
Nt—i = 2Nt—i. Considering the initial condition Nq = 2, 
we obtain N t = 2 t+1 and n v (t) = 2*. Since U t is a tree, we 
have the following relation Et = Nt — 1 = 2 t+1 — 1. On the 
other hand, at arbitrary step t > 1, the addition of each 
new node leads to only new edge, thus n e (t) — n v (t) = 2 t 
for all t > 1. 



3 Topological properties 

Topological features of a network are of fundamental sig- 
nificance to understanding the complex dynamics taking 
place on it. Here we focus on four important characteris- 
tics of the tree Ut, i.e., degree distribution, average path 
length, betweenness distribution, and degree correlations. 

3.1 Degree distribution 

The degree is the simplest and most intensively studied 
characteristic of an individual node. The degree of a node i 
is the number of edges in the whole network connected to i. 
The degree distribution P(k) is defined as the probability 
that a randomly selected node has exactly k edges. Let 
ki(t) denote the degree of node i at step t. If node i is 



added to the network at step U, then by construction, 
ki{ti) = 1. In each of the subsequent time steps, a new 
node will be created connected to i. Thus the degree ki{t) 
of node i satisfies the relation 

h(t) = ki(t - 1) + 1. (1) 

Considering the initial condition kiiti) = 1, we obtain 

ki(t) =t-ti + l. (2) 

Since the degree of each node has been obtained explicitly 
as in Eq. ([2]), we can get the degree distribution via its 
cumulative distribution [5] 

oo 

P cum {k) = ]T P(k'), (3) 
fe'=fc 

which is the probability that the degree is greater than 
or equal to k. An important advantage of the cumula- 
tive distribution is that it can reduce the noise in the 
tail of probability distribution. Moreover, for some net- 
works whose degree distributions have exponential tails: 
P(k) ~ e~ fc / K , the cumulative distribution also has an 
exponential expression with the same exponent: 

oo oc 
Pcum(k) = ~ E e ~ k ' /K ~ e ~ k/K - W 

This makes exponential distributions particularly easy to 
detect experimentally, by plotting the corresponding cu- 
mulative distributions on semilogarithmic scales. 

Using Eq. ©, we have P cum {k) = £~ =fc P(fc') 
= P (*' < t = t + 1 - k). Hence 

Pcum(k) = = 2 *+i " = 2 ^ ^ ( 5 ) 

which decays exponentially with k. Thus the determinis- 
tic uniform recursive tree is an exponential network, and 
has a similar form of degree distribution as its stochastic 
version — the random uniform recursive tree 1251 . 



3.2 Average path length 

Average path length (APL) means the minimum number 
of edges connecting a pair of nodes, averaged over all node 
pairs. It is defined to be: 

Q 

dt = N t (N t - l)/2> (6) 

where St denotes the sum of the total distances between 
two nodes over all pairs, that is 

S t = E d i ,3 ' ( 7 ) 

i¥=3 

where di.j is the shortest distance between node i and j. 
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Fig. 1. Illustration of the deterministic uniform recursive tree, showing the first five steps of growth process. 



Let /?ncw an d ^oid represent the set of nodes created 
at step t or earlier, respectively. Then one can write the 
sum over all shortest paths St in network Ut as 

S t = E dl >J> 

' (8) 
where the third term is exactly St-i, i.e., 



In the infinite network size limit (t — > oo), 



(14) 



(9) 



which means that the average path length shows a log- 
arithmic scaling with the size of the network, indicating 
a similar small-world behavior as the URT [30 and the 
Watts-Strogatz (WS) model gg]. 



By construction, we can obtain the following relations for 
the first and second terms: 

]T d i , j =Nl 1 + 2S t - 1 , (10) 

tki = + N t -i(Nt-i ~ I)- (11) 

Substituting Eqs. Q, dTU]) and (HI]) into Eq. (© and con- 
sidering N t — 2 t+1 , the total distance is obtained to be 

S t = 4S t _ 1 + 2A/ t 2 _ 1 -A/ t _ 1 

t-i t—i 
= 4* So + 2 ^ 4 t ~ 1 - l iVf - ^ 4*- 1 - i 7V i 



i=0 

At I 



i=0 



= 2< • 4 f + 2 J 

Inserting Eq. (fl"2]) into Eq. ([6]), we have 

- _ 2(2i ■ 4* + 2*) _ <-2 t+1 + l 
* _ 2*+ 1 (2*+ 1 - 1) ~~ 2*+i - 1 ' 



(12) 



(13) 



3.3 Betweenness distribution 



Betweenness of a node is the accumulated fraction of the 
total number of shortest paths going through the given 
node over all node pairs [491150], More precisely, the be- 
tweenness of a node i is 



E 



CTjk 



(15) 



where ajk is the total number of shortest path between 
node j and k, and <Jjk{i) is the number of shortest path 
running through node i. 

Since the considered network here is a tree, for each 
pair of nodes there is a unique shortest path between 
them [51, 52,53, 54J . Thus the betweenness of a node is sim- 
ply given by the number of distinct shortest paths passing 
through the node. Then at time t, the betweenness of a 
r-generation-old node v, which is created at step t — T + 1, 
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denoted as bt(r) becomes 



b t (r) = 0\ [N t - (OJ + 1)] + 



E 



6>™ + 1 
2 



(16) 

where OJ denotes the total number of descendants of node 
v at time t, where the descendants of a node are its chil- 
dren, its childrens children, and so on. Note that the de- 
scendants of node v exclude v itself. The first term in 
Eq. (fT6|) counts shortest paths from descendants of v to 
other vertices. The second term accounts for the shortest 
paths between descendants of v. The third term describes 
the shortest paths between descendants of v that do not 
pass through v. 

To find bt(r), it is necessary to explicitly determine the 
descendants Q\ of node v, which is related to that of v's 
children via [5511541 



oi = y: (o{ + 1) 



Using Q\ = 0, we can solve Eq. (|T7|) inductively, 
01 = 2 T_1 - 1. 



(17) 



(18) 



Substituting the result of Eq. ([HJ and N t = 2 t+1 into 
Eq. (fT6|) . we have 



Mr) 



-)t + T 



■)t + l 



2 2(r-l) + f 5 



(19) 



which is approximately equal to 2 t+T for large t. Then the 
cumulative betweenness distribution is 

which shows that the betweenness distribution exhibits a 
power law behavior with exponent 75 = 2, the same scal- 
ing has been also obtained for the URT [53] and the m = 1 
case of the Barabasi-Albert (BA) model [55] describing 
a random scale-free treelike network [5T][52] . Therefore, 
power-law betweenness distribution is not an exclusive 
property of scale- free networks. 



3.4 Degree correlations 

An interesting quantity related to degree correlations 56 
is the average degree of the nearest neighbors for nodes 
with degree k, denoted as k nn (k) [57ll58ll59] . When k nn (k) 
increases with k, it means that nodes have a tendency 
to connect to nodes with a similar or larger degree. In 
this case the network is defined as assortative .60,61] . In 
contrast, if k nn (k) is decreasing with k, which implies that 
nodes of large degree are likely to have near neighbors with 
small degree, then the network is said to be disassortative. 
If correlations are absent, k nn (k) = const. 

For the deterministic uniform recursive tree, we can 
exactly calculate k nn (k). Except for the initial two nodes 



generated at step 0, no nodes born at the same step, which 
have the same degree, will be linked to each other. All links 
to nodes with larger degree are made at the creation step, 
and then links to nodes with smaller degree are made at 
each subsequent steps. This results in the expression for 
k = t + l-U (U > 1) 



k n n{k) — 



1 



n v (ti)k{ti,t) 
t!-=t 



tj=ti-i 

E 

t;=o 



ra„(i-)/c(4t) 



n v (ti)k(t'i,t) 



=u+i 



(21) 



where k(ti,t) represents the degree of a node at step t, 
which was generated at step U. Here the first sum on the 
right-hand side accounts for the links made to nodes with 
larger degree (i.e. t ■ < U) when the node was generated at 
t i . The second sum describes the links made to the current 
smallest degree nodes at each step t\ > ti. 

After some algebraic manipulations, Eq. (I2ip is sim- 
plified to 



knn{k) 



t t% t — (— 3 t% 



t + l-ti t + l-U 



(22) 



Writing Eq. ([2"2"| in terms of k, it is straightforward to 
obtain 



k+l 



k ■ 2* 



(23) 



Thus we have obtained the degree correlations for those 
nodes born at U > 1. For the initial two nodes, each has 
a degree of k = t + 1, and it is easy to obtain 



k nn (k 



t + 



1 



(24) 



t'=o 



From Eqs. (I23p and ([24]) . it is obvious that for large net- 
work (i.e., t — ► 00), k nn (k) is approximately a linear func- 
tion of k, which shows that the network is assortative. 



4 Eigenvalues and eigenvectors of the 
Laplacian matrix 

As known from section 2, there are 2 t+1 vertices in Ut- we 
denote by Vt the vertex set of Ut, i.e., Vt — {vi, Kg, ... , v 2 t+i}. 
Let At = [ay] be the adjacency matrix of network Ut, 



where a, 



a 7 j 



1 if nodes i and j are connected, 



aji = otherwise, then the degree of vertex is defined as 
dv i = J2jev t a ii ■ Let Dt = diag(d t , 1 ,d V2 ,.. . , <4 2t+1 ) repre- 
sent the diagonal degree matrix of U t , then the Laplacian 
matrix of Ut is defined by L t — D t — A t . For an arbitrary 
graph, it is generally difficult to determine all eigenvalues 
and the corresponding eigenvectors of it Laplacian ma- 
trix, but below we will show that for Ut one can settle 
this problem. 
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4.1 eigenvalues 

We first study the Laplacian spectra of U t making use of 
an iterative method [62] • By construction, it is easy to find 
that the adjacency matrix A t and diagonal degree matrix 
D t obeys the following relations: 



and 



D, 



Af_i If_i 
It i 



D t _i+I t _i 

It_i I ' 



(25) 



(26) 



where each block is a 2* x 2* matrix and I t _i is identity 
matrix. Thus, according to the definition of Laplacian ma- 
trix, we have the recursive relation between L t and L t _i 
as 



D t - A t = 



D t _i+I t _i- A t -i h 
-lt-i It 

Lt-i + It— i — It-i 
-It-i It-i 



(27) 



Then, the characteristic polynomial of L t is 
g t (x) = det(xl t - L t ) 

V it-i (x-i)i t _i 

/ _ 1 L_\ T. . 

= det 




where the elementary column operations of matrix have 
been used. According to the results in [53], we have 

9t{x) 



det y\x - 1 - — - J It-i -Lt-i) ■ det((x - l)I t _i ) 

It-l-Lt.i). (29) 



l) 2 * -det( (x- 1 1 



x- 1 



Thus, gt(x) can be written recursively as follows: 



g t (x) = (x- l) 2 ■ g t -i(f(x)), 



(30) 



where f(x) — x — 1 — This recursive relation given 
by Eq. (|30[) is very important, from which we will deter- 
mine the complete Laplacian eigenvalues of Ut and their 
corresponding eigenvectors. Notice that gt-i(x) is a monic 
polynomial of degree 2*, then the coefficient of l/(x— l) 2 
in g t -i{f{x)) is 1, and hence 1 is the constant term of 
gt{x). Consequently, 1 is never an eigenvalue of L t . 

Note that Ut has 2* +1 Laplacian eigenvalues, and all 
these eigenvalues are distinct, which will be shown below. 
Let these 2 t+1 eigenvalues are A^, A|, . . . , A 2t+1 , respec- 
tively. For convenience, we presume A* < A| < • ■ ■ < AL+i 



and denote by A t the set of these eigenvalues of Ut, i.e., 
A t —{X\, A 2 , . . . , A 2t+1 }. 

From Eq. (|30p . we have that for an arbitrary element 



in A t -i, say A* 1 € At-i, both solutions of x — 1- 



A' 1 are in At 



are m Zlj. In fact, equation x — 1 — — ^t- — A* 1 is 
equivalent to 



(A 4 *- A +2)x + Aj- l = 0. (31) 

We use notations A* and A* +2 t to represent the two solu- 
tions of Eq. ([3"T]) , since they provide a natural increasing- 
order of the Laplacian eigenvalues of Ut , which can be seen 
from below argument. Solving this quadratic equation, 
its roots are obtained to be A* = ri(A- _1 ) and \\ +2 t = 
r 2 (A* _1 ), where the function r*i(A) and r 2 (A) satisfy 



r-i(A) = i(A + 2- v / A 2 +4), 



r 2 (A) = -(A + 2 



V\ 2 



(32) 
(33) 



Substituting each Laplacian eigenvalue of Ut-i into Eqs. 
([32]) and (f33|) , we can obtain the set A t of Laplacian eigen- 
values of U t . Since A — {0, 2}, by recursively applying the 
functions provided by Eqs. ([3"2"]l and ([3^]) . the Laplacian 
eigenvalues of Ut can be determined completely. 

It is obvious that both r\ (A) and r 2 (A) are monotonously 
increasing functions. On the other hand, since the in- 
dependent variables A here are greater than or equal to 
zero, so ri(A) > and r 2 (A) > 0. Furthermore, fi(A) = 
(A + 2-VA 2 + 4)/2 = (A-VA 2 + 4)/2 + l, where the term 
(A — vA 2 +4)/2 is clearly less than zero, thus r\[X) < 1. 
Similarly, we can show that r 2 (A) > 2. Thus, for arbi- 
trary fixed A', ri(A) < r 2 (A') holds for all A. Then we 
have the following conclusion: If the Laplacian eigenvalues 
set of U t -i is A t -i = {A* -1 , A 2 -1 , . . . , A^ 1 }, then solving 
Eq. (|31|) one can obtain the set of the Laplacian eigen- 
values of U t to be A t = {X\,X\, ■ ■ ■ , A| t , A* t+1 , . . . ,A 2t+1 } 
where < A* < A 2 < . . . < A 2t < 1 < 2 < A 2t+1 < 
A 2t+2 < ... < A 2t+ i. Therefore, all the 2' +1 Laplacian 
eigenvalues of Ut are different, which has never been pre- 
viously reported in other network models thus may have 
some far-reaching consequences. 



4.2 eigenvectors 

Similarly to the eigenvalues, the eigenvectors of Lt follow 
directly from those of L*_i. Assume that A is an arbitrary 
Laplacian eigenvalue of U±, whose corresponding eigenvec- 
tors is v, then we can solve equation (A If — Lt)u = to 
find the eigenvector v. This equation can be also written 
as 

/ \ t x \ ( (A - l)I t -i - Lf_ x I t _i A ( vi\ 
(^-^>={ I t _ x (A-l)I i _ 1 J( vW2 J= ' 

(34) 

where V\ and v 2 are two components of v. Eq. ([34]) leads 
to the two following equations: 

((A-l)I t _i-L t _i)w 1 +i;a=0, (35) 
V! + (A - l)v 2 = 0. (36) 
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Resolve Eq. ([36J to find 
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l 



V2 



-Vl 



A- 1 

Substituting Eq. ([37]) into Eq. ([35]) we have 
1 



A- 1 



A - 1 



t-i — M-i 



=0, 



(37) 



(38) 



which indicates that Vi is the solution of Eq. (J3SJ) while 
V2 is uniquely decided by V\ via Eq. ([37)) . 

From Eq. ((30|) in preceding subsection, it is clear that if 
A is an eigenvalue of Laplacian matrix L t , then / (A) = A — 
1 — must be one eigenvalue of L t _i . (Recall that if A = 

A| 6 A, then /(A*) - A, 4 " 1 for i < 2 4 , or /(A*) = A>: 2 \ 
for i > 2*). Thus, Eq. (J2HD together with Eq. J30) shows 
that t>i is an eigenvector of matrix L t _i corresponding to 
the eigenvalue A — 1 — determined by A, while 



v = 



l 

A-l 



(39) 



is an eigenvector of Lt corresponding to the eigenvalue A. 

Since for the initial graph Uq, its Laplacian matrix Lo 
has two eigenvalues and 2 with respective eigenvectors 
(1,1) T and (1,— 1) T . By recursively applying the above 
process, we can obtain all the eigenvectors of the deter- 
ministic uniform recursive tree lit- 



5 Conclusion and discussion 

We have investigated a deterministic model for the uni- 
form recursive tree, which is constructed in a recursive 
way. The model is actually a deterministic variant of the 
intensively studied random uniform recursive tree. We have 
presented an exhaustive analysis of many properties of the 
considered model, and obtained the analytic solutions for 
most of the topological features, including degree distribu- 
tions, average path length, betweenness distribution, and 
degree correlations. Aside from their deterministic struc- 
tures, the obtained statistical characteristics are equiva- 
lent with its corresponding random URT. Consequently, 
the DURT may provide useful insight to the practices as 
URT. 

Also, we have performed a detailed analysis of the com- 
plete Laplacian spectra and eigenvectors of DURT by us- 
ing the methods of linear algebra and graph theory. We 
have fully characterized the spectral properties and eigen- 
vectors for DURT. We have shown that all the eigenval- 
ues and eigenvectors of the Laplacian for DURT can be 
directly determined from those for the initial graph. In- 
terestingly, all the Laplacian eigenvalues have only one 
multiplicity. To the best of knowledge, this property has 
never been reported for other previously studied model. 
Finally, it should be mentioned that although we have ex- 
amined only a specific model, the methods used here can 
be applicable to a larger type of networks. 
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